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The iteration procedure of supersymmetric transformations for the two- 
dimensional Schrodinger operator is implemented by means of the matrix 
form of factorization in terms of matrix 2x2 supercharges. Two different 

> 

00 ■ types of iterations are investigated in detail. The particular case of diagonal 



initial Hamiltonian is considered, and the existence of solutions is demon- 



m 

O ■ strated. Explicit examples illustrate the construction. 

o 

C^' 1. Introduction 

Supersymmetric Quantum Mechanics (SUSY QM) [1], [2] has become a well-known pow- 
erful tool for the investigation of problems of modern non-relativistic Quantum Mechanics. 
The main achievements of SUSY QM during the last two decades involved isospectral pairs 
of one-dimensional systems. In particular, the conventional SUSY QM (with supercharges 
of first order in derivatives) was generalized [3], [4] to Higher order SUSY QM (HSUSY QM) 
with supercharges n— th order polynomials in momentum. This generalization induces the 
deformation of SUSY QM algebra: the anticommutator of supercharges is now a polynomial 
in the Superhamiltonian. The class of pairs of isospectral systems which are connected (in- 
tertwined) by components of supercharge is thus enlarged [4]. Two classes of second order 
supercharges were found: reducible, which can be written as a chain of consecutive first 
order transformations, and irreducible second order transformations. In both cases the 
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intertwined Hamiltonians are expressed in terms of only one arbitrary function [4]. For a 
(admittedly nonexhaustive) list of references devoted to one- dimensional HSUSY QM see, 
in particular, [5]. 

It is also useful to recall the equivalent terminology related to the classic mathematical 
paper by G.Darboux [6]. The conventional first order SUSY transformations coincide with 
the so-called Darboux transformations, which are very useful in Mathematical Physics (see, 
for example, [7] for nonlinear evolution equations). Then n— th order HSUSY transforma- 
tions correspond to the well-known Crum-Krein transformations [8], where the eigenfunctions 
of n-th iterated Hamiltonian are obtained in terms of n given eigenfunctions of the initial 
Hamiltonian by means of "determinant expressions". The second iterations of SUSY trans- 
formations for matrix one-dimensional potentials were considered in [9] , recently the general 
form of matrix Crum-Krein transformations was constructed [10]. 

Much less attention was paid in the literature to SUSY QM of two-dimensional 
systems. Here we list some partial achievements which have been obtained in this field. 
The direct multi-dimensional generalization of the standard one-dimensional SUSY QM was 
elaborated in [11], [12] for arbitrary space dimensions. In the simplest two-dimensional case 
the supersymmetrization of a given (scalar) Schrodinger Hamiltonian involves not only a 
second scalar Hamiltonian but also a matrix 2x2 Schrodinger operator. The spectra of 
these three components of the Superhamiltonian are interrelated, and their wave functions 
are connected by the components of supercharges. This approach was used successfully 
[13], [14] to find the spectrum and the wave functions of Pauli operators describing spin 
1/2 particles (with arbitrary values of the gyromagnetic ratio) in a broad class of non- 
homogeneous external electromagnetic fields in terms of spectra and eigenfunctions of two 
scalar Hamiltonians. 

In the series of papers [15] direct intertwining relations between a pair of scalar two- 
dimensional Hamiltonians were studied. A wide variety of particular solutions was found, 
leading to new integrable quantum systems which are not amenable to separation of vari- 
ables. This approach allowed to introduce [16] two new methods for the study of the two- 
dimensional Schrodinger equation: SUSY— separation of variables and shape invariance. 
Multi-dimensional SUSY QM for an arbitrary oriented Riemannian manifolds was studied 
in [17]. In a very recent paper [18] the Hamilton- Jacobi approach of Classical Mechanics was 
applied to study classical two-dimensional supersymmetric models and their quantization. 
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The present paper addresses the problem of the construction of a chain of consecutive 
super symmetric (generalized Darboux) transformations for two-dimensional Schrodinger sys- 
tems allowing for a deeper insight into isospectrality in multidimensional cases. Let us remark 
that the only early attempt to construct such iterations failed (see Section 4 of the first paper 
in [11]). Even the very possibility of making iterations for two-dimensional systems remained 
a task to be solved. 

The new matrix supersymmetric factorization of the two-dimensional Pauli operator, 
recently proposed in [14], can provide a new chance to construct iterations of first order 
(reducible) supersymmetric transformations. Indeed, in contrast to the method of quasifac- 
torization [11], where the vector components of supercharges include scalar superpotential, 
the matrix factorization of [14] involves the components of supercharges with matrix super- 
potentials. Therefore all elements of the SUSY intertwining relations are now represented 
by 2 x 2 matrix differential operators. 

The organization of the paper is the following. In Section 2 the chain of of SUSY trans- 
formations is performed by gluing the corresponding components of two successive Super- 
hamiltonians. The generalization of this prescription includes an additional global unitary 
transformation (Section 3). In Section 4 we analyze the particularly interesting case when the 
initial Hamiltonian is diagonal, and one supersymmetrizes a scalar Hamiltonian by build- 
ing partner Hamiltonians via second order SUSY transformations. This construction, in 
reverse order, can be also used for diagonalization of a class of two-dimensional 2x2 matrix 
Schrodinger operators or for matrix systems of coupled second order differential equations. 

2. Second iteration by gluing 

Before introducing the new form of SUSY factorization for two-dimensional systems, 
proposed in [14], let us start from the main formulas for the conventional [2] one- dimensional 
SUSY QM with partner Hamiltonians and their eigenfunctions: 



#(°) = q + q 



d 2 + V^(x); H^\x) 



(1) 





The operators q ± are defined as 




d + W{x); 



q- = ( q +y = + d + W(x); 



d = d/dx, 
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and they intertwine the components (1), (2) of the Superhamiltonian: 

if<V = g+ifW; q~H^ = H^q~. (3) 
Then, up to normalization factors, 

*SHx) = <r*i 0) (*); *i 0) (*) = q + *SHx). (4) 

The intertwinings (3) are the most important relations of the SUSY QM algebra: 

{Q+, Q-} = H; (Q + ) 2 = (Q-) 2 = 0; [H, Q±] = 0, (5) 
where the Superhamiltonian H and supercharges are: 

h= \ H(0) ° h Q + =(Q-y=\ " : i. (6) 




o 

This procedure of supersymmetric (Darboux) transformations from H^°\ ^f^\x), E n to 
H^ 1 \^ n 1 \x),E n was iterated [4] by reducible supercharges of second order in derivatives. 
This iteration leads to a polynomial deformation of SUSY algebra. 

A two-dimensional generalization of standard SUSY QM [11] organized (in the Super- 
hamiltonian) one matrix 2x2 and two scalar Hamiltonians 6 . Each scalar Hamiltonian 
is separately intertwined with the matrix Hamiltonian, but the spectra of the two scalar 
Hamiltonians are not related. This construction was based on the quasi-factorization of all 
three components of the Superhamiltonian in terms of the first order differential operators 
qf = =F<9/ + (dix( x )) and-^ pf = e/fcgj; l,k — 1,2. A similar quasi-factorization with new 
supercharges, used in [11] to build an iteration of SUSY transformations, turned out to be 
unsuccessful since several necessary constraints were too complicated to be implemented [11]. 
In addition, it led to a growth of the matrix dimensionality of the new Hamiltonians. 

In the present paper we use the method [14] based on the matrix factorization of com- 
ponents of the two-dimensional Hermitean Superhamiltonian 9 : 

if(°) = q+q- = -A (2) + {Wi) 2 + a 1 {d 2 W l -d l W 2 )-az{d 1 W 1 + d 2 W 2 ) (7) 



Tor the case of an arbitrary dimensionality d of coordinate space the Superhamiltonian contains (d+ 1) 
components H^ n \ n = 0, 1, d with matrix dimensionality C%, i.e. dimensionality of the space with definite 
fermionic occupation number n (see details in [12]). 

•^The summation over repeated indices (i = 1,2) is implied here and below. 

9 Slightly different notations are chosen here in comparison with [14]. 
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ifM = q - q + = -^) + (Wif + a l (d 2 W 1 + d 1 W 2 ) + a z (d 1 W l -d 2 W 2 ) (8) 
9,= ^- (i = l,2) A = c\ 2 + d 2 2 

The 2x2 components g r± of 4x4 supercharge ^ (see Eq.(6)) are matrix differential operators 
of first order in derivatives h : 

q ± = Tdi - icr 2 d 2 + a x W 2 + a 3 W 1 q + = (g") f , (9) 

where W\ )2 are real, and a,i {i = 1,2,3) are standard Pauli matrices. By construction, these 
operators intertwine the Hamiltonians (7), (8) according to the standard relations (3) of 
SUSY algebra. The relations (4) - (6) of one-dimensional SUSY QM are satisfied as well, 
with the proviso that we deal now with 2x2 matrix differential operators q ± , H^°\ 
and two-component wave functions Isospectrality of the (matrix) Hamiltonians 

H(°\ H^> holds, as usual, except for the possible zero modes of q ± . 

To iterate supersymmetric transformations we impose the ladder equation (see [3], [4]) 
thus gluing the lower component of the Superhamiltonian H with the upper component 
if(°) of the next Superhamiltonian H, defined by (7), (8) with W replaced by W : 

if W = #(°) + C, C = const{real). (10) 

Then the second order supercharge Q + = q + q + gives the following intertwining relations: 

tf(°)Q+ = Q+(hW + C). 

The conditions (10) can be rewritten as a nonlinear system of differential equations: 

Wf = W? + C; (11) 

d 2 w x + d ± w 2 = d 2 w ± - d 1 w 2] (12) 

d 1 W 1 -d 2 W 2 = -d{W x -d 2 W 2 . (13) 

It may be convenient to express these ladder equations in terms of complex function of two 
mutually conjugated complex variables 

W(z, z*) = Wi(z, z*) + iW 2 (z, z*); z = x x +ix 2 z* = x x - ix 2 , (14) 



''An analogous form of matrix supercharges was used also in [19] in the context of ParaSUSY QM. 
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as follows: 



\W(z,z*)\ 2 = \W(z,z*)\ 2 + C; 
dW(z,z*) = -dW(z,z*); 
d=— = -(81-182); 8 



(15) 
(16) 
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(d 1 +id 2 ). 



dz 2 V 1 zn dz* 2' 

The general solution of Eq.(16) can be written by means of an arbitrary complex function: 

W(z,z*) = -8F(z,z*); W(z,z*) = 8F(z,z*); F(z, z*) = F ± (z, z*) + iF 2 (z, z*) (17) 

with real F\ and F 2 . 

In terms of Fi(z, z*), F 2 (z, z*) Eq.(15) reads: 

(ftFiX&Fa) - (8^(8^) = C. (18) 

Its general solution is a general solution of the homogeneous equation summed to a par- 
ticular solution of the inhomogeneous one. One can check that the general solution of the 
homogeneous equation requires that F 2 (x±, x 2 ) depends on its arguments via Fx(xi, x 2 ) only: 



F 2 (x u x 2 ) = $ 



Fi(xi,x 2 ) 



F(z, z*) = x 2 ) + iF 2 (x u x 2 ) 



(19) 



(and, of course, vice versa). In order to find the particular solution of inhomogeneous 
equation (18) one has to introduce new variables yi,y 2 such that y 2 = x 2 and yi coincides 
with one of the solutions (19), e.g. y\ = Fi(x\,x 2 ). In new variables Eq.(18) becomes: 

dF 1 {x 1 ,x 2 ) dF 2 (y u y 2 ) 



dxi 

and its solution for F 2 in terms of F\ is: 



dy 2 



F 2 (x 1 ,x 2 ) = C J dy 2 



' dF 1 (xi,x 2 ) 
dx 1 



where the partial derivative is taken with x 2 = const, and after integration over y 2 the initial 
variables x±,x 2 are to be reinserted. So the general solution of (18) is the sum: 

'9F 1 (xi,x 2 )\ _1 



F 2 (x 1: x 2 ) = $ 



F 1 (x 1 ,x 2 ) 



+ Cj dy 2 \ 



dxi 



(20) 



Thus two matrix two-dimensional Hamiltonians (components of new Superhamiltonian) 
(see expression (7)) and (H^ + C), which can be written as: 



ifW + C = # (0) + 2a 1 8 1 8 2 F 1 (x 1 ,x 2 ) - 2a 3 d 1 d 2 F 2 (x 1 ,x 2 ), 



are intertwined by the differential matrix operator (component of new supercharges): 

Q+ = q +q+ = dl-dZ + 2ia 2 d 1 d 2 + 
1 



+ 4 
+ °i 



(d 2 F 2 ) 2 + (d 2 Fi) - (d.F.Y - (d 1 F 1 ) 
(8^)8, - (d 2 F l )d l - 1(3? + d 2 2 )F 2 



+ 



- r 2 



(9 1 F 1 )(9 2 F 1 ) + (9 1 F 2 )(9 2 F 2 ) 



+ 



1 



(9 2 F 2 )9! - (d.F^d, + ^(d 2 + d 2 2 )FA 



+ 0"3 

The SUSY algebra has now a polynomial form: 

{Q+, Q-} = - C); (Q + ) 2 = (Q-) 2 = 0; [if, Q± 



0. 



This algebra differs essentially from the deformed SUSY algebra for a class of two-dimensional 
systems investigated in a series of papers [15] in so far as it does not incorporate any non- 
trivial symmetry operators (no central extension): thus we have no positive indications of 
integrability of the systems H^°\ H^. 

Higher order iterations. In order to perform the next iteration of supersymmetric 

~(o) 

(Darboux) transformations one has to glue IF > with the first component H of a new 
Superhamiltonian. This procedure follows the previous one: superpotentials are expressed 
in terms of derivatives of the same arbitrary complex function F(z, z*) : 

W = W x +iW 2 = +dF(z, z*) = -OF; (21) 
W = W l +iW 2 = +dF(z,z*). (22) 

The last equality in (21) implies that in terms of an arbitrary function f(z,z*) : 

F(z, z*) = -df(z, z*)- t F(z, z*) = +df(z, z*). 

It is difficult to find the general solution of the nonlinear equations (15) and the analogous 
one for W, W, but it is possible to find particular solutions with certain simplifying ansatzes. 
For example, explicit solutions can be obtained for C = C = and f(z,z*) = |0(^)| 2 . 



3. Iterations by gluing and global transformations 

We now modify the procedure (namely, Eq.(10)) of the previous Section by gluing the 
Hamiltonians and up to a global unitary transfomation U. Such a procedure 
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allows to preserve the main details of the scheme, but it might extend the class of intertwined 
Hamiltonians. 

Without loss of generality, we may restrict ourselves* to U — icr 3 : 

if W = UH^U- 1 + C, U = ia 3 . (23) 

In this case (15), (16) must be replaced by: 

\W(z,z*)\ 2 = \W(z,z*)\ 2 + C; (24) 
d 1 (W 2 -W 2 ) = ^(Wi + ^i); (25) 
d 1 (W 1 + W 1 ) = d 2 (W 2 -W 2 ). (26) 

The Eqs.(25), (26) can be solved explicitly: W and W are interrelated by an arbitrary 
holomorphic function G(z) : 

W(z, z*) + W*(z, z*) = G(z); G(z) = d+ iG 2 ; d 1 G 1 = d 2 G 2 ; d ± G 2 = -d^. (27) 

The remaining nonlinear equation (24) for G can be written as: 

G\ + G\- 2W 2 G 2 - 2W 1 G 1 + C = 0. (28) 

It can be solved explicitly as a linear algebraic equation for superpotentials W\,W 2 in 
terms of an arbitrary holomorphic function G(z). For an arbitrary G(z) one can choose 
also an arbitrary function W 2 (z,z*) and find the functions Wi(z,z*), Wi :2 (z,z*) from Eqs. 
(28), (27). Thus one obtains a class of Hamiltonians (7) for which the second order matrix 
supersymmetric transformations can be performed. The resulting Hamiltonian 

ff (!) + c = ff {0) + 2^ • [d x W 2 - d 2 W x + d 2 G^j + 2a 3 • d 1 G l (29) 

is intertwined with the initial Hamiltonian by the operator 

Q+ = q +a 3 q + , 

and is therefore isospectral (up to zero modes of Q~^) to H^°\ 

Example with hidden symmetry. Among many possible examples we consider 
the simple peculiar case when the resulting ff ^ coincides with the initial up to a 

constant: 

H (°) = H (1) + C. (30) 
*From now on we will refer to such global unitary transformation as "rotation" even disregarding i. 
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This can be achieved by imposing (see (29)) G(z) = iaz and 

1 C 1 C 

W 2 (z, z*) = -(ap H ) cos0 + 6(p) sin0; WAz, z*) = -(ap H ) sin0 + 0(p) cos0, 

2 p 2 p 

where p, are polar coordinates, a is a constant and O(p) is an arbitrary function. The 
Hamiltonian (30) then becomes: 

H (o) = _ A (2) + hap+ -f + (6(p)) 2 - aa, - -(pe(p))V 3 . (31) 

P P 

We stress that though the potential in (31) depends on p only, its matrix structure in general 
prevents a standard separation of variables, typical for scalar problems. From the intertwin- 
ing (now commutation) relations [H^ , Q^] = one can find two mutually commuting 
Hermitean symmetry operators: 

R~ = ^(Q + - Q-) = tadf, (32) 

1 1 C 

R + = -(Q + + Q~) = -a 3 H {0) + iaa 2 pd p + a 3 -ap{ap + -) - a ia pQ{p), (33) 

with column-eigenfunctions, common to H^°\R ± : 

The first one, R~ , reflects the obvious symmetry of under rotations., but the second 
symmetry operator (33) realizes certain "hidden" symmetry and has the property: 

(R + ) 2 + (R-) 2 = H(°\H<® - C). 

One can notice that the spectrum of the Hamiltonian involves double m <-> — m degeneracy, 
and in addition a possible double degeneracy (for each value of m) associated with a column- 
function ifj m (p). A pedagogical example illustrating the considerations above is obtained from 
(31) with 0(p) =0 and C = : it is exactly solvable after suitable rotation and leads to a 
pair of decoupled radial oscillators with frequiency a and a relative ground state energy shift 
2a. 

Higher order iterations. In general, the procedure of gluing with a global rotation 
enlarges essentially the class of Hamiltonians for which our method is applicable. We now 
demonstrate this considering two options to construct third order transformations: 

I. The first option reads: 

~ ~ ~(o) „ ~(i) 

# (0) — > //(i) = #(0) + c — > + C = a 3 H a 3 + C + C — >H +C + C, 



where — ► denotes the appropriate first order matrix SUSY transformation induced by 
q + ,q + ,q + , and = represents the gluing. The first gluing is of the kind described in Sec- 
tion 2 (Eq.(lO)), and the superpotentials W(z,z*) and W(z,z*) satisfy Eqs.(17), (18). The 
second gluing is described in Eq.(23), and W(z,z*), W(z,z*) satisfy Eqs.(27), (28). 
Eq.(28) can be rewritten as: 

(G(z)dF* + G*(z*)dF) - \G\ 2 -(7 = 0, (34) 

where F(z,z*) still satisfies (18). We now restrict ourselves to the case C — 0, for which 
Eq.(18) has trivial particular solutions with purely real or purely imaginary F. Then, e.g. 
for real F, Eq.(34) can be solved by replacing z, z* by: 

^ I 2G{z) + I 2G^) ] T ^ l {l ' 2G{z)~ I 2G^])' (35) 
since Eq.(34) in new variables reads: 

d t F ± = \G\ 2 + C. 

A solution for F^. 

F 1 = T(t) + J \G\ 2 dt + Ct 

depends on two arbitrary functions - one holomorphic (G(z)) and one real (T(r)). The 

F = F 1 so obtained generates W and W when inserted into (17) and (27) (in terms of the 

~(i) 

initial variables xi, x 2 ). Thus one obtains a class of Hamiltonians H (°\ (H +C + C), which 
are intertwined (and are isospectral) by the third order operator Q + = q + q + a 3 q + . 

2. The second option includes two gluings, both with a 3 -rotations: 

~ ~ ~(o) „ ~(i) 

H ^ — > H (1) = a 3 H (0) a 3 + C — > H {1) + C = a 3 H a 3 + C + C — > H +C + C. 

The superpotentials W(z, z*) and W(z, z*) has to be found from equations of the form 
(28), and then W(z,z*) and W(z, z*) from equations of the form (27). One can obtain the 
general solution of these equations purely algebraically in terms of two arbitrary holomorphic 
functions G(z), G(z) : 

, G 2 G - \G\ 2 G - CG - CG* , x 

W(z, z*) = U = , 36 

GG - G*G* 

W(z, z*), W(z, z*) being derived, again algebraically, from two equations of form (27). Thus 
one obtains a new class of Hamiltonians H (0 >, (H + C + C), which are intertwined by the 
third order Q + = q + a 3 q + a 3 q + , and are therefore isospectral up to zero modes of Q ± . 
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4. Second order supersymmetrization of scalar two- 
dimensional Hamiltonians. 



In this concluding Section we will discuss particular cases of Sections 2 and 3, when the 
initial Hamiltonian is constrained to be diagonal. The condition of diagonality of (7) 
reads: 

d 2 W 1 -d 1 W 2 ^0 & W 2 = d 2X (x 1 ,x 2 ); W 1 = d lX (x 1 ,x 2 ), (37) 

where x( x ii x 2) is an arbitrary real function. In terms of this function one has = 
diag(H$, H { 22 ] ) with 

H[? = -A< 2 ) + (d lX ) 2 - (fl? + dl) X ; H% = + 2(9? + dl) X , (38) 

which is intertwined (and therefore, isospectral, up to zero modes of supercharges (9)) with 
non-diagonal Hamiltonian The latter is obtained by inserting (37) into (8). From now 
on we study the iteration of SUSY transformations for the initial Hamiltonian (38), i.e. the 
compatibility of conditions (37) with the iteration algorithm. 

Simple gluing. For the case C = in (10) the combination of Eq.(17) with Eq.(37) 
can be solved in terms of an arbitrary real function £(xi,x 2 ) : 

F 1 (x 1 ,x 2 ) = {dl-dD^x^ (39) 
F 2 { Xl ,x 2 ) = -29 1 9 2 e(x 1 ,x 2 ); (40) 
X(xi,x 2 ) = (d 2 1 +dMxi,X2)- (41) 
The last condition to be taken into account is the interrelation (19) between Fi and F 2 : 

(d 2 2 - %)S{x u x 2 ) = ffr&tixuxz)]. (42) 

Of course, one can not find the general solution of this nonlinear functional-differential 
equation, but one can find a variety of particular solutions. The simplest choice / = const 
in (42) is trivial, since the corresponding Hamiltonian if ^ can be directly diagonalized by 
rotation. Some less trivial solutions of (39), (40), (42) can be found by requiring that the 
l.h.s. of Eqs.(39), (40) are e.g. F { = \£ k ; Fi = expA^, etc. 

For the case C^O one can try similar techniques as used for the solution of Eq.(18) to 
solve now the equation (18) and (37) together with (17). The resulting equation reads: 

(d 2 1 -d 2 2 )F 2 = 2d 1 d 2 F 1 . 
11 



Particular solutions of this equation can be found by suitable ansatzes like, for example: 

F 1 (x 1 , x 2 ) = x ig (x 2 ) + f 2 {x 2 ); $(Fi) = aF? + bF 1 + c+^ r , 

where g(x 2 ), f(x 2 ) can be determined by solvability conditions, the functional $(iq) has been 
introduced in (20), and a, b, c, d are constants. One can also enlarge the class of solutions by 
suitable nonsingular linear transformations of Xi,x 2 . 

Gluing with rotation. In this case the nonlinear equation (28) takes the form: 

(d lX )G 1 + (d 2X )G 2 = l/2[G\ + G\ + C], 

or, in terms of variables z, z* : 

{G{z)d + G?{z*)&) X {z,z*) = \G(z)\ 2 + \c. (43) 

It has the form similar to Eq.(34) of the previous Section except for the fact that x( z , z *) is 
a real function. Therefore, the general solution of (43) is expressed in terms of t, r of (35) 
via an arbitrary (real) function A(r) : 

X(z,z*) =A(r) + / \G{z)\ 2 dt+ l -Ct. (44) 

The initial diagonal Hamiltonian H^> from (7) and (37) reads: 

H (o) = _ A (2) + 4(dx)(dx) -4{ddx) -<r 3 = 

- "3-(G' + G*'+|^), 

while the final can be expressed as: 

+ C = H (0) - 2a 1 9iG 2 + 2a 3 d 1 G 1 . 

As a particular example one can consider the case of G = z 1 ! 2 which leads to a confining 
singular Hamiltonian with a singularity 1/p and growing asymptotically as p with a nontrivial 
azimuthal dependence. 

Finally, we note that conversely the intertwining between diagonal and non-diagonal 
by second order matrix supercharges can also be used to diagonalize a class of matrix 
Schrodinger operators (or system of differential equations of second order with non-diagonal 
matrix coupling), in analogy with similar procedures already used in [13], [14], [20]. 
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